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Beam-recoil transferred polarizations for the exclusive ep — » e' K~^ti, E° reactions have been mea- 
sured using the CLAS spectrometer at Jefferson Laboratory. New measurements have been com- 
pleted at beam energies of 4.261 and 5.754 GeV that span a range of momentum transfer from 
0.7 to 5.4 GeV^, invariant energy W from 1.6 to 2.6 GeV, and the full center-of-mass angular range 
of the K'^ meson. These new data add to the existing CLAS K^ A measurements at 2.567 GeV, 
and provide the first-ever data for the K'^TP channel in electroproduction. Comparisons of the data 
with several theoretical models are used to study the sensitivity to s-channel resonance contributions 
and the underlying reaction mechanism. Interpretations within two semi-classical partonic models 
are made to probe the underlying reaction mechanism and the ss quark-pair creation dynamics. 
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PACS numbers: 13.88.+e, 14.40.aq, 14.20.Gk, 14.20.Jn 



I. INTRODUCTION 

An important requirement to better understand the structure of the nucleon is to map out its spectrum of excited 
states. However, deciphering the data to understand the resonance excitations has been hmited both by the data itself 
and the current state of existing theories. Ideally we should expect the fundamental theory of the strong interaction, 
quantum chromodynamics (QCD), to provide a prediction for the nucleon excitation spectrum. However, due to the 
non-perturbative nature of QCD at the relevant energies, this idea has not yet been fully realized. Thus we have looked 
instead to effective models of QCD, such as constituent quark models, to gain some insight. Present quark model 
calculations of the nucleon spectrum have predicted more states than have been seen experimentally 1[. This has 
been termed the "missing" resonance problem, and the existence of these states is tied in directly with the underlying 
degrees of freedom of the nucleon that govern hadronic production at moderate energies 2| ■ 

Most of our current understanding of nucleon resonances comes from reactions involving pions in the initial and/or 
final states. Koniuk and Isgur suggested that the missing states might be revealed in decays to channels where mesons 
other than pions or multiple pions are in the final state [3| . Indeed, there are indications from theory that some missing 
states have a similarprobability of decaying into channels such as loN , r]N, ttttN, and KY {Y = A, S) compared to 
the nN channel Q, 0]. As baryon resonances have large widths and are often overlapping, studies of different final 
states provide important complementary cross checks in quantitatively understanding the contributing amplitudes. 

In this work we study the electroproduction of strange final states. While electromagnetic production of KY final 
states has a much lower cross section than hadronic production reactions, the use of an electromagnetic probe has 
a distinct advantage, namely that all electromagnetic quantities in the reaction amplitude can be straightforwardly 
expressed in the context of quantum electrodynamics. Furthermore, in addition to the different coupling constants 
compared to the ttN channel (e.g. qkny vs. g^NN), the study of the exclusive production of KY final states has 
another advantage in the search for missing resonances. The higher masses of the kaon and hyperons, compared to 
their non-strange counterparts, kincmatically favor a two-body decay mode for states with masses near 2 GeV. Not 
only is this situation advantageous from an experimental viewpoint, but this also happens to be the mass region where 
the majority of the missing resonance states are expected to exist 

Although the two ground-state hyperons have the same valence quark structure (uds), they differ in isospin, such 
that intermediate N* resonances can decay strongly to KA final states, while both N* and A* decays can couple to 



KY^ final states. Existing studies of N* KA, JCS and A* — > ifS decays have not yet provided extensive or precise 
information on the N*,A* KY couplings. To date, the Particle Data Group (PDG) only lists four N* states 
with known couplings to KA and no N* states are listed that couple to KI] [5,]; only a single A* state is listed with 
coupling strength to if E. The current landscape as given by the PDG for N*, A* — > KY is given in Table [J 
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TABLE I: PDG listings for the coupling of TV* (A*) states below 2 GeV to KA and KT. (KT.) [g. The Rating column gives 
the PDG star rating for the TV* states and B.R. indicates the branching ratio. 



Theoretically, there has been considerable effort during the past two decades to develop models for KY photo- and 
electroproduction. However, the present state of understanding is limited by a lack of precision data (Ref. 6| contains 
a brief review). Model fits to the cross section data are generally obtained at the expense of many free parameters, 
which makes it difficult to provide precise constraints. Moreover, cross section data alone are not sufhcient to fully 
understand the reaction mechanism, as they represent only a portion of the full amplitude response. In this regard, 
measurements of spin observables are essential for continued theoretical development in this field. Fits to the limited 
available data lead to ambiguities and model dependence in interpreting the results. Polarization data can provide 
for improved constraints on the model parameters, increasing their discriminatory power and allow for a quantitative 
measure of whether or not new resonance states are required to explain these and other hyperon production data. 
One main issue involves discriminating resonant states from the non-resonant background and from effects caused by 
final-state interactions or channel- couplings instead of TV* and A* contributions [if. 

CLAS at Jefferson Laboratory (JLab) has provided photoproduction K^A and if+S" recoil polarization data from 
the proton [8]. In addition, beam-recoil polarization transfer data from CLAS have been published for both K^ A and 
_ftr+S" photoproduction 9] and K^A electroproduction 10| reactions on the proton. Data such as these that span 



both a wide energy and angular range and are essential to disentangle the resonant and non-resonant contributions 



to the KY spectrum 7 



llj. This has been demonstrated in several recent ampHtude-level analyses with channel 



couplings based on photoproduction data 



12, 



13 



14 



15 



16 



17 



181 ligj . Further progress is expected as data 



with broad coverage and smaller experimental uncertainties are made available (which includes new CLAS data with 



linearly polarized photon beams and polarized targets 



21 



In this work, we focus on measurements of spin transfer from a longitudinally polarized electron beam to the 
ground-state hyperons produced in the reactions p{e,e'K^)A and p{e, e' K^)'S'~' at beam energies Ei, of 4.261 and 
5.754 GeV. This work represents a higher-statistics follow-up to the first data presented by CLAS in the A'+A channel 
for an electron beam energy of 2.567 GeV 10|, where the transferred A polarization was studied as a function of the 
invariant energy W and cos6'^™' (the center-of-mass angle). The transferred polarization data for the A'+E" 
final state included here represent the first-ever published data for this observable in electroproduction. 

From the polarization data in Ref . [IQ] , the ratio of the longitudinal to transverse structure functions (Tl /ct^ for the 

A'+A final state at 6*^"'- = 0° was extracted for several W points near 1.8 GeV and ^0.7 GeV'^ 23]. These results 

indicated a ratio that was systematically smaller than previously published results using a Rosenbluth separation 

I i 

performed in Hall C at Jefferson Laboratory [2J|, albeit with large statistical uncertainties. In fact, the data were 
consistent with zero within the experimental uncertainties, which would imply a small longitudinal structure function, 
and hence, a small longitudinal coupling of the virtual photon. The results of Ref. j23| are expanded upon in this work 
with larger data sets that reduce uncertainties in the extrapolation to f?^™- ~ 0°. The new data presented include 
three data points near 11^=1.9 GeV with an average of ~1.6 GeV^ and three data points near W—2.0 GeV with 
an average of ~2.5 GeV^. 

Using a semi-classical partonic framework, the CLAS polarization data in Ref. [10[ were shown to support a 



description where the spin properties of the quark-pair creation operator might be responsible for the observed trends 
in the A polarization. This framework indicated that the quark-pair creation operator dominating the reaction 
produces the ss pair with spins anti-aligned. This finding, if confirmed, has important implications since many, if not 



most, calculations of hadron spectroscopy use a ^Pq operator to calculate the transition to the final-state particles 251. 
In this work, the angular distribution of the transferred polarization is studied with greater precision than in Ref. 10 1 
and compared against two semi-classical partonic models that lead to quite different predictions regarding the A'+A 
reaction mechanism and the quark-pair creation dynamics. The first is the model from Ref. 10|. The second assumes 



the reaction proceeds from an ss quark-pair with the quark spins aligned. The main differences between the models 
are discussed and a possible experiment to discriminate between them is proposed. 
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The organization for the remainder of this paper is as follows. In Sections HIl and Hill the theoretical models to be 
compared with the measurements are briefly introduced and the relevant formalism for the polarization measurements 
is provided. Section [TVl gives a detailed description of how the polarization is extracted and Section IVl gives details 
regarding the analysis cuts and corrections to the data. Section IVII details the sources of systematic uncertainty. 
Section fVIII contains the physics results, with the presentation of the if+A and K^YP polarization transfer data in 
Sections IVII Al and IVII Bl respectively, the new a^/aT extraction in Section lVlI CI and comparisons of the data to the 
newly developed partonic models in Section IVII Dl Finally, we present a summary of this work and our conclusions 
in Section Ivml 

II. THEORETICAL MODELS 

While the QCD description of quark interactions and pair creation is well accepted at high energies, the situation 
is considerably more complex in the low-energy nucleon resonance region due to the non-perturbative nature of the 
theory. In order to arrive at any theoretical expectations for the transferred polarization, effective models must be 
employed that ultimately represent approximations to QCD. This analysis highlights three different theoretical model 
approaches. The first is a traditional hadrodynamic model, the second is based on kaon Regge trajectory exchange, 
and the third is a hybrid Regge plus resonance approach. 

A. Hadrodynamic models 

Hadrodynamic models provide a description of the reaction based on an effective Lagrangian constructed from 
tree-level Born and extended Born terms in the s, t, and u reaction channels (see Fig. [1]). The Born diagrams include 
the exchange of the proton, kaon, and ground-state hyperons, while the extended Born diagrams include the exchange 
of the associated excited states. This description of the interaction, which involves only first-order terms, is sensible 
as the incident and outgoing electrons interact rather weakly with the hadrons. A complete description of the physics 
processes requires taking into account all possible channels that could couple to the initial and final states, but the 
advantages of the tree- level approach are to limit complexity and to identify the dominant trends. The drawback 
in this class of models is the large number of hadrons that can contribute in the intermediate state of the reaction. 
Depending on which set of resonances a given model includes, very different conclusions about the strengths of the 
contributing diagrams may be reached. 
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(a) 

s-channel 



(b) 
t-channel 




FIG. 1: Tree-level diagrams contributing to the KY reactions: (a) s-channel exchanges, (b) t-channel exchanges, and (c) 
u-channel exchanges. 



The hadrodynamic model employed in this work is from Mart and Bennhold 26| (referred to here as MB). The 
s-channcl terms included in this model are listed in Table [TTl The coupling strengths have been determined by a 
simultaneous fit to low-energy K^p —> 'jY and 7*^*^^ —> K^Y data, by adding the non-resonant Born terms with 
a number of resonances, leaving the coupling constants as free parameters. The coupling constants are required to 
respect the limits imposed by SU(3), allowing for a symmetry breaking at the level of about 20%. In this model, the 
inclusion of hadronic form factors leads to a breaking of gauge invariance that is restored by the inclusion of counter 



terms [26!]. The model has been compared to the CLAS photoproduction [8 



27[ and electroproduction data [6] and 



provides for a fair description of those results, although no CLAS data were employed in the model fits. 
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TABLE II: List of s-channel resonant terms included in the MB model [26[ and the Regge plus resonance (RPR) model 0] 
included in this work. Note that the RPR model has two variants that include either a _Di3(1900) or Pii(1900) state. 



For AT+A production, the MB model includes four baryon resonance terms. Near threshold, the steep rise of the 
cross section is accounted for with a core set of N* states: 5*11(1650), Pn (1710), Pi3(1720). To explain the broad 



bump in the energy dependence of the cross section seen by SAPHIR 



281 and CLAS [61,18 



27|, the MB model includes 



a spin-3/2 £)i3(1900) resonance that was predicted in the quark model of Capstick and Roberts 



l| to have a strong 



8 



coupling to the fC+A channel, but which was not well established from existing pion-production data. For i^+S° 
production, the MB model includes the core N* states and the A* resonances 5*31(1900) and P3i(1910). The model 
also includes if* (892) and 7^1(1270) exchanges for both KY final states, but does not include any u-channel diagrams. 

The N* states S'ii(1650), Fii(1710), and Pi3(1720) are the only states listed by the Particle Data Group jsl with 
coupling strengths to KA (see Table |l| . While the relevance of these core states in the K^A reaction has 

long been considered a well-established fact, this set of states falls short of reproducing the experimental results in the 
region below 14^=2.0 GeV. Furthermore, two recent analyses have called the importance of the Pii(1710) state into 



question 



16 



Beyond the core states, the PDG lists a two-star Pi3(1900) state as the sole established N* near 
1900 MeV. However, with a 500-MeV width, it appears unlikely that this state by itself can explain the structure(s) 
visible in the CLAS and SAPHIR cross sections, unless its parameters are significantly different than those given by 
the PDG. This has led to suggestions of a new (unconfirmed) N* state in this mass region (e.g. the Z3i3(1900) state 
in the MB model). However, the analysis of Saghai 



29( 1 ■ using the same data sets employed for the MB model fits, 
concluded that by tuning the u-channel background processes involved in the K^A reaction, the need to include any 
states beyond the core set was removed. Note that the investigation of contributing N* states to the KY reactions 
has typically been limited to spin j < 5/2 due to the expectations that higher-spin resonances do not significantly 
contribute to the reaction dynamics. 



16 



30| 



Moving beyond tree-level approaches to consider recent multipole and coupled-channels models has not led to 
dramatic new insights to the N* spectrum. The multipole analysis by Mart and Sulaksono 
coupled-channels models of Julia-Diaz et al. [17] and Sarantsev et al. 18| (which all employ CLAS photoproduction 



3l|, as well as the 



data in their fits), claim that a _Di3(1900) state is required by both the CLAS and SAPHIR 773 — ^ K^A data. However, 



the coupled-channels model of Ireland et al. 



points to a Pii(1840) state as a more likely candidate (although one 



or more of Sn, Pu, P13, D13 are not ruled out). The fits of Julia-Diaz et al. 17] suggests a third 5ii resonance might 



also be playing a role, while Sarantsev et 



18| also require (in addition to a Di3(1900)) the presence of a Pii(1840) 



18| by Avisovich 



and another D13 state at 2170 MeV. An extension of the coupled-channels model of Sarantsev et al. 
et al. [19| , which was the first model to include the CLAS photoproduction hyperon polarization transfer observables 
Cx and Cz 9|, concluded that a Pi3(1900) state was also required to satisfactorily fit the data. 

In the recent fits of the jp K^I]'^ data, all N* resonances found to be necessary to fit the K'^A data have been 
included. However, the existing X+E'^ database is much smaller than the K'^A database, with significantly larger 
statistical uncertainties. Even with this situation, the recent coupled-channels models 



17 



18 



19[ indicate important 



9 

resonant contributions to the K+YP final state from the N* states Pii(1840), Di3(1870), Pi3(1885), and Di3(2170), 
and from the A* states ^35(1905), ^33(1940), and ^37(1950). 

Each different model has ambiguities that can be better constrained only by incorporating better quality data 
or including new experimental observables. Comparison of the models to the data can be used to provide indirect 
support for the existence of the different baryonic resonances and their branching ratios into the strange channels, as 
well as improved constraints on the phenomenology of the different strangeness production reactions. 



B. Regge and regge plus resonance models 



Our KY electroproduction data are also compared to the Reggeon-exchange model from Guidal, Laget, and Van- 



derhaeghen [33| (referred to here as GLV). This calculation includes no baryon resonance terms at all. Instead, it is 
based only on gauge-invariant t-channel K and K* Regge-trajectory exchange. It therefore provides a complementary 
basis for studying the underlying dynamics of strangeness production. It is important to note that the Regge approach 
has far fewer parameters compared to the hadrodynamic models. These include the K and K* form factors (assumed 
to be of a monopole form) and the coupling constants qkyn and qk'YN (taken from photoproduction studies). The 
GLV model was fit to higher-energy photoproduction data where kaon exchanges dominate and extrapolated down to 
JLab energies. Furthermore, the use of Regge propagators eliminates the need to introduce strong form factors in the 
background terms, thus avoiding the gauge-invariance issues associated with traditional effective Lagrangian models. 
The GLV Regge model reason ably accounts for the strength in the CLAS A'+A differential cross sections and 



separated structure functions [a, |27|. Although the reasonable performance of a pure Regge description in this 
channel suggests a ^-channel dominated process, there are obvious discrepancies with the data, indicative of s-channel 
strength. In the K^YP channel, the same Regge description significantly underpredicts the differential cross sections 



and separated structure functions g, |27[. The fact that the Regge model fares poorly when compared to the K^YP 
data is indicative that this process has a much larger s-channel content compared to K production. 

The final model included in this work is based on a tree-level effective field model for A and Y'^ photoproduction from 
the proton. It differs from traditional isobar approaches in its description of the non-resonant diagrams, which involve 
the exchange of K and K* Regge trajectories. A selection of s-channel resonances are then added to this background. 
This "Regge plus resonance" model (referred to here as RPR) tI has the advantage that the background diagrams 
contain only a few parameters that are constrained by high-energy data where the t-channel processes dominate. In 
addition to the kaonic trajectories, the RPR model includes the s-channel resonances 511(1650), Pii(1710), Fi3(1720), 



and Pi3(1900) (see Table HI]) . The model also includes either a D 



In detailed comparisons with the separated structure functions y, |3J] and transferred polarization data from CLAS 
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1^1900) or Pii(1900) state in the K+K channel. 



10|, only the £'i3(1900) assumption could be reconciled with the data, whereas the Pii(1900) option could clearly 



be rejected 



In the K+YP channel, four A* states, D33(1700), S'3i(1900), P3i(1910), and P3i(1920), have been 



included (see Table HI)) . 



III. POLARIZATION FORMALISM 



A. Polarization component definitions 



The differential cross section for kaon electroproduction can be written as the product of a virtual photon flux 
factor r^, and the kaon virtual differential cross section, expressed in the kaon center-of-mass (cm.) frame as 



da 



da,, 



dne'dn^^'dEe' dn'j;^^- 



(1) 



The most general form for the differential cross section of a kaon from a proton target, allowing for a polarized electron 
beam, target proton, and recoil hyperon, is given by 35 1 



day 

dn^ 



R 



pa 



eRl" + c+ (Xt cos $ + Xt sin $) 



pa 



eC'R^'^ cos 2$ + 'R^° sin 2$) 



hc^ {"R^ cos $ + Xt' sin $) + /icoi?^^' ■ 



Pa 



(2) 



The R^'^ terms represent the response functions that account for the structure of the hadronic system and, in general, 
are functions of Q^, W , and cos^^;-™- only. The superscripts a and (5 refer to the target and hyperon polarization 
axes, respectively, and the c and s superscripts indicate a cosine or sine dependence on the angle $, where $ is the 
angle between the electron and hadron planes. Here e is the transverse polarization of the virtual photon, h is the 
electron-beam helicity, and K, is the ratio of the momentum of the kaon to the virtual photon in the cm. frame The 



factors c± are given by -\/e(l ± e) and co = \/l — e^. Fig. [2] defines the angles of the scattering process in the cm. 
system. 

The operators Sa and Sp project out the target polarization vector in the {x,y,z) system and the hyperon polar- 
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electron scattering 
plane 




FIG. 2: Kinematics for KY electroproduction defining tiie cm. angles and coordinate systems used in the analysis. This figure 
shows a positive <E> angle. 

ization in the (x', y' , z') system, respectively (see Fig. [2]). The (x, y, z) system is defined such that z is along the three 
momentum transfer q direction and y is normal to the electron-scattering plane. The {x' , y' , z') system is defined such 
that z' is along the kaon momentum vector and y' is normal to the hadronic plane. 

In the case where there is no beam, target, or recoil polarization (/i, a, (3 — 0), Eq.([2]) reduces to 



ao=IC [R^ + eR"j^ + c+R% cos $ + eRfj, cos 2$] . 



(3) 



For the case of a polarized-electron beam incident on an unpolarized target producing a polarized recoiling hyperon, 



Eq.([2) becomes 



36| 



dcr. 



"— = ao(l + hALT' + Px'x' ■ S^. + Py,y' ■ Sy, + P,,z' ■ S,,) 



(4) 



where Alt' = ^c^Rfj., sin <f> is the polarized beam asymmetry defined in terms of the fifth response function R^t' ■ 
Each of the recoil- hyperon polarization components can be split into a beam-helicity-independent part , called 
the recoiZ polarization, and a beam-helicity-dependent part P/, called the frans/errerf polarization. The components 
of the hyperon polarization vector can be written as Pi = Pf + hP'^. The three recoil polarization components are 
given in terms of the response functions in the (x', y' , z') system as 



P°, 



p3 



— (c+Rl^ sin $ -t- e R^^ sin 2$ 
do V 

/C 



R 



yo 



eR 



v'o 



c+RH'S. cos $ + eRl^ cos 2$ 



— \c+Rlt sin $ + ei?TT sin 2$ 
(Jo V 



(5) 
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and the three transferred polarization components are written in the (a;', y' , z') coordinate system as 



P'^, = — (c-Rl^, cos $ + coR^TT' 
an V 



CTO 

P' = — c_i?^'" sin $ 



P',, = — (c^Rl^, cos $ + . (6) 

Co V / 



To accommodate finite bin sizes and to improve statistics, our analysis sums over all $ angles. The ^-integrated 
polarization components (represented by the V symbol) in the {x\ y' , z') system are given in Table Hill In performing 
the (f> integration, the polarization components P", , T'^, , and V'y, are equal to zero. In Table [1111 the term Kj ~ 
l/(i??P + ei?oo). 



Coordinate System 



' x' 





V , 

' x' 


KiCqR^°, 


■po 

y' 


Ki{R%° + eRy^°) 


V , 

v' 





pO 





V, 

' z' 


KjCqR^^/ 


(x, y, z) Coordinate System 


-pO 

' X 





V 


\Kic^{Ri^, cos6»^™- - + R{^, sin 6*^" ) 


pO 


ii<r/c+(i?|'^cos6l^™- + i?iTsin6l^™) 


V 





pO 





v. 


KiCo{-R^^, sin6'^"- + iJf,'^, cos 6*^" ) 



TABLE III: Polarization observables integrated over $ in the two coordinate systems used in this work. 



To define the polarization observables in the {x, y, z) coordinate system, shown in Fig. [2l the components defined 
for the {x',y',z') system in Eqs.® and ((6]) must undergo a simple transformation that involves a rotation of 
about y' , followed by a rotation of $ about z'. The ^-integrated recoil and transferred polarization components in 
the {x,y,z) system are defined in Table Hill 



B. (Tl/ (Tt ratio 



The A polarization transfer data can be used to extract the ratio of the longitudinal-to-transverse structure functions 
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33] and for cos O'jf^- ^ 1 [6| 



10| taken at a beam 



~ oi^joT at cos^^™- = 1. This ratio has been previously measured at cos6'|^'"' = 1 
using the Rosenbluth separation technique. Our previously published polarization transfer results 
energy of 2.567 GeV were used to extract R^r at cos6'^™- = 1 for three points of W and Q|. In a similar way, R„ 
can be extracted from the 4.261 and 5.754 GeV CLAS data. 
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In parallel or anti-parallel kinematics (cos^^"*- = ±1), the z' and z components of the transferred polarization 
integrated over $ (given in Table UlTl) reduce to 



V'. ^±V'^± 



± 



(7) 



where the plus (minus) sign is associated with the parallel (anti-parallel) kinematics case and cr^ = (7^ + 6(7^. 

The response functions used to express the ^-integrated components of V',, and V'^ (see Tabl ejnl)) can be written 



in terms of the Chew, Goldberger, Low, and Nambu (CGLN) amplitudes 



3811 as shown in Ref. 



For the case of 



^c^m. — 0°, it can be shown that R^^, ~ R^, and Eq.([7]) can be rewritten as 



vl, = ri = 



(8) 



Inverting this form and rearranging, the ratio Ra at cos0|^™' = ±1 can be written as 



Ra 



Oh 
(Tt 



1 



(9) 



While the V'^, and V'^ data presented here do not include data points at ^l^™- — 0°, an extrapolation to cos6'^™- = ±1 
can be performed as shown in Section IVII CI 



IV. POLARIZATION EXTRACTION 



A. Decay angular distributions 



The A decays weakly into a pion and a nuclcon with the decay nucleon constrained to move preferentially in the 
direction of the hyperon spin. In the A rest frame, the decay nucleon angular distribution is given by jsQ ] 



dN 



COS0^^), 



(10) 



where Pa is the A polarization and is the angle between the polarization axis and the decay-nucleon momentum 
in the A rest frame. In this work we focus solely on the A pn^ decay (B.R.=64%) and explicitly replace Opf^ with 
9^^ ■ The A weak decay asymmetry parameter Q!a has been measured to be 0.642±0.013 5|. 
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The S° decays into a 7 and a A (branching ratio 100%). A TP with polarization Ps will yield a decay A that 
retains some of the polarization of its parent. As shown in Ref. 40|, Pa = — -jPs on average for the decay A in its 
rest frame. For the case of a final-state the A rest frame can be calculated only if four particles are detected in 
the final state. In addition to the detection of the electron, kaon, and decay proton, either the decay pion of the A 
or the decay 7 from the S*^ must be detected. Due to the small CLAS acceptance for a four particle final state, only 
three final-state particles were detected. In Ref. [9] it has been shown that the polarization of the daughter A from 
the decay can be measured without boosting the detected proton to the reference frame of the A. The value of 
the effective weak decay asymmetry parameter was determined to be as = —0.164, or in terms of the A weak decay 
constant, as — — 0.256aA, thus reduced from the value — 0.333aA. This value is independent of kinematics. 

As the electron beam is not 100% polarized, the hclicity term h in the hyperon polarization must be replaced by 
the average longitudinal electron-beam polarization Pf, as 

PY=P^+PbP{-. (11) 

Combining the expressions from Eqs. pop and pip , the decay proton angular distribution for the two different beam 
helicity states can be written 

'^^^^N^\l + o.{P^±PMcose^% (12) 



d cos „p 



where q;a = 0.642 for the A analysis and as = —0.164 for the YP analysis. 



B. Asymmetry approach 

The transferred hyperon polarization was extracted using the acceptance-corrected yield asymmetry for the two 
different electron beam helicity states of the form 

. N+ -N- 

This asymmetry is formed from the cos 9^^ yields for the three spin-quantization axes of the decaying hyperon (in 
either of the coordinate systems defined in Fig.[2|). The terms and represent the acceptance-corrected decay 
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proton yields in a given kinematic bin. The helicity-gated yields are given by 



iV±(cos0^^) = E = aoE[l ± P^Alt' + a{P^ ± PtP^) cos^p«^], 



(14) 



where £ represents the CLAS detection efficiency, which is assumed to be a helicity-independent function, and includes 
the CLAS acceptance function and the beam-target luminosity factors. 

As discussed in Section IIIIl this analysis was performed by integrating the decay proton yields over all $ angles to 
maximize the statistical precision of the measurement. In this case the measured yield asymmetry becomes 

_ " ao£[l + PbALT' + a{P^ + PbPy) cos - ao£[l - PtALT' + a{P^ - P^Pi,) cos 

/o'" ao£[l + PbALT' + a(po + P^P^) cos 0^^^]d$ + " cj^£[l - P^Alt' + a{P^ - PbP^) cos eRF\d<^ 

After some simplification, this expression can be written as 

_ /o " ao [PbALT' + aPbPy cos e^^]d<^ 

— ^TZ ■ ( J^t) ) 

/o"ao[l + aPO cosf?^f^]d$ 

As the Alt' term is proportional to sin^/co, it integrates to zero for all choices of spin-quantization axes. Similarly, 
the term containing Py integrates to zero along the {x,x') and (z, z') axes using the definitions in Section [IIII 

Considering our coordinate system choices (see Fig. [5]) , it turns out that the ^-integrated asymmetries can only 
be non-zero along the {x^x') and {z,z') axes. This can be seen as J^^ aoPyd^ ~ along both the {y,y') axes, 
given the polarization definitions in Section [IIII In other words, when performing the $ integration, the polarization 
components PP ,,n, V?,^ and V'/, are all constrained to be zero. Thus the only possible non-zero asymmetries 
for our coordinate systems will be j,/) and A(^^ ,j.fy Along these special axes, the ^-integrated asymmetries can be 
written as 



A = aPb 



CaoPi-d^' 



cos^^^. (17) 



The quantity in brackets is equivalent to the V polarization terms in Table IIIIl and represents the ^-integrated 
hyperon transferred polarization, where a separate asymmetry is computed for each spin-quantization axis. We 
therefore extract the non-zero transferred hyperon polarizations with respect to the different quantization axes from 
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the asymmetries by fitting 



= aAPbV'f, cos 9^^ or As = asP^T'^ cos 0^^ . (18) 



In forming the asymmetry of Eq. (jl3p . the decay proton helicity-gated yields are sorted for each kinematic bin of 
interest. We used an event-by-event weighting factor to correct the yields for the detector acceptance. The asymmetry 
method used in this analysis is relatively insensitive to the detailed form of the CLAS acceptance function. This 
discussion is contained in Section [VD I 

C. Hyperon polarization and statistical uncertainty 

For the general case where a given hyperon sample Y is contaminated by particle misidentification events and events 
from the tail of the hyperon Y' , the helicity asymmetry can be written in terms of its individual contributions as 

_ (jV+ + N + + N+,) {Ny + Ny, + N,-,) 



Ny + Ny + Nbck 



+ 



''bck ^''bck . ^bck 



Ny ~^ Ny 

where Ny, Ny', and N^ck refer to the number of counts from the hyperon of interest, the tail of the other hyperon, 
and from the background (mostly pions misidentified as kaons), respectively, within the Y identification cuts (see 
Section IVE)) . If we define Fy = Ny /Ny and F^ck — Nbck/Ny, we can write 

, _ Ay +AyFy + AbckFbck ,„„s 

' l + Fy+Fb.k • ^ ^ 

It was observed in this analysis that the pion background asymmetry under both hyperon peaks is consistent with 
zero, thus the term associated with Ajjck in Eq.((20|) is set to zero (see Section IV Fp . The "pure" asymmetries, Aa and 
As, are given by Eq. (jl8|) . thus the measured asymmetry can be written in terms of the transferred polarizations, as 
well as the measured polarization within the hyperon mass cuts 'Pmeas 

A ^yP^ncosor + -y'P^'Py'Cose^'Fy 

1 + ry + tbck 
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where cihPbP'jj^f,g^g is the slope extracted from the fit of the cosO^^ distribution with respect to a given spin- 
quantization axis. 

We found that the contamination of E° hyperons within the A identification cuts is consistent with zero (see 
Section IV Ep . so for Y — A and F' = E", Eg. ipTjl can be rearranged to get 

K=K.eas{^+nck). (22) 

For y = S° and Y' = A, Eq.([2T|l can be rearranged to get 

= — [-PLasi^ +Fa+ Ftck) - n^A] . (23) 

as 

Performing standard error propagation, the statistical uncertainty for V'^ is 

= [(1 + Pbck){SVUeas? + iPLasfSFL] , (24) 

where the individual uncertainties are given by 



Here (5(slope) is the uncertainty in the slope from the fit of the cosOp^ distribution. We find that the dominant 
contribution to the A polarization uncertainty is due to the uncertainty in T^^eas- 
Similarly, the statistical uncertainty for is 

SV^ = ^ [(1 + Fa + n,,r{6V'^,,^r + {V'meas " Kf^Fl + {VUeas5F,ck? + [F^SV'^f] , (26) 

as 

where the individual uncertainties are given by 



= + (27) 



SFbck = Fbck\l ^ + and (28) 

V ^^bck JVS 

- (29) 
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The dominant contribution to the E'' statistical uncertainty arises due to the uncertainty in V^^^g. All other terms 
are at least a factor of 4 smaller in size. 



D. Depolarization factor 

Our formalism defines the polarization transfer as the ratio of the hyperon polarization to that of the electron 
beam. However, the electron interacts with the hadronic system through the exchange of a virtual photon. Thus 
the true "beam" polarization is given by the product Pi,D{y), where D{y) accounts for the polarization loss from the 
incident beam electron to the virtual photon. There are a number of ways to express the factor D{y). One form is 
given by 41 1 



where y = E^* / Ei, is the relative energy transfer to the target proton and = ctl/ct- 

With this accounting, the hyperon polarization can be rewritten from Eq.((ni) as Py = + PbD{y)P^. This 
would lead to slightly modified forms of the asymmetries in Eq. (jl8p with the hyperon polarizations scaled by a 
factor of 1/D{y). This re-expression of the hyperon transferred polarization allows for a more direct comparison for 
experiments performed at different beam energies. Perhaps, more importantly, it allows for a more direct comparison 
of electroproduction and photoproduction data sets where the "trivial" depolarization factor is accounted for in the 
electroproduction data. 

Having made this distinction in the possible convention choice for the hyperon transferred polarization, we have 
decided not to account for it in this work, following instead the procedures in Ref. lj|. The main reason for this 
choice is to avoid introducing a model-dependent uncertainty into our quoted polarizations. Our studies have shown 
that with different hadrodynamic models, the variation in D{y) due to variations in can be up to 20%. 



42|. 



Shown in Fig. [3] are predictions of the depolarization factor for several of our kinematic bins using the MB model 
It is seen that in our kinematics {D(i/)) ^ 0.6. The polarization transfer from the virtual photon to the hyperon is 
therefore 67% larger on average than that for the beam electron to the hyperon. When considering the polarization 
data in Section rVIIl one must take into account the depolarization factor when comparing to other data and to theory. 
The theory calculations shown in Section IVIII match the data, namely they show the product of the depolarization 
factor and the polarization. 
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Q" (GeV ) 

FIG. 3: Calculations of the depolarization factor at 5.754 GeV using the MB model [13] in representative data sorts (a) 

summing over all and cos6'J>'"' and (b) summing over all W and cos^f^""'. 



V. DATA ANALYSIS 



A. The CLAS detector 



43|. 



All of the data shown in this analysis were collected using the CLAS spectrometer located in Hall B at JLab 
The main magnetic field of CLAS is provided by six superconducting coils, which produce an approximately toroidal 
field in the azimuthal direction around the beam axis. The gaps between the cryostats are instrumented with six 



identical detector packages, as shown in Fig. [4] Each sector consists of three sets of drift chamber (DC) packages 



44| 



to determine the trajectories of the charged particles, Cherenkov counters (CC) ,_45|| for electron identification, scin- 



tillator counters (SC) 



46| for charged particle identification, and electromagnetic calorimeters (EC) 



47| for electron 



identification and detection of neutral particles. A 5-cm long liquid-hydrogen target was located in the center of the 
detector on the electron beam axis. 

To reduce the electromagnetic background resulting from M0ller scattering off atomic electrons in the target and 
the target cell, a small normal-conducting toroidal magnet (called the mini-torus) was placed symmetrically about 
the target inside of the first DC package. This magnetic field sweeps M0ller electrons out of the detector volume. 
A totally absorbing Faraday cup, located at the end of the beam line, was used to determine the integrated beam 
charge passing through the target. The efficiency of detection and reconstruction for stable charged particles in the 
fiducial regions of CLAS is greater than 95%. The solid angle coverage of CLAS is approximately Stt sr. The polar 
angle coverage for electrons ranges from 8° to 45°, while for hadrons it is from 8° to 142°, with an angular resolution 
of 66, (50 ^ 1 of better than 2 mr. The CLAS detector was designed to track particles having momenta greater than 
roughly 200 MeV with a resolution Sp/p in the range of 0.5 to 1%. 
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FIG. 4: (Color online) Three dimensional view of the CLAS detector with the different subsystems labeled. A single sector of 
the detector has been cut away to enable a view of the inner subsystems. The diameter of the CLAS detector is ~5 m and it 
is ~8 m long. 

The large acceptance of CLAS enabled us to detect the final-state electron and kaon, as well as the proton from 
the decay of the A hyperon. Hyperon identification with CLAS relies on missing-mass reconstructions of the reaction 
e + p e' + + X. In this section, details are provided on our procedures for particle identification, the cuts used 
to isolate the K'^ K and if+S" final states, the hyperon spectrum fitting procedures, and other cuts and corrections. 



B. Data set information 

The data were taken with typical electron beam currents of 5 nA at a luminosity of 1 x lO'^^ cm^^s^^. The 
CLAS event readout was triggered by a coincidence between a Cherenkov counter and a calorimeter detector in a 
single sector, generating an event rate of ~2 kHz. The main CLAS torus had its polarity set such that negatively 
charged particles were bent toward the electron beam line. The electron beam was longitudinally polarized, with the 
polarization determined by a coincidence M0ller polarimeter. Beam polarization measurements were taken at regular 
intervals throughout the running periods and measured a stable electron beam polarization of 70%. 

The data in this paper were collected as part of the CLAS running periods elc in early 1999 and el-6 in late 
2001/early 2002. The elc running period included data with beam energies of 2.567 GeV (previously published in 
Ref. 



10[) and 4.261 GeV acquired at several different field settings of the main CLAS torus. 4.261 GeV represents 
the luminosity-averaged beam energy for data taken with electron beam energies of 4.056, 4.247, and 4.462 GeV. 
Combining the data sets is justified given the relatively small spread in the virtual photon polarization parameter 
e among the different energies. The el-6 running period was taken with an electron beam energy of 5.754 GeV. 
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Information regarding the different run periods, including the total number of triggers, the approximate W and 
ranges of the data, the number of hyperons in the different analyses detected through the e'K^p final state, and the 
average beam polarization, is contained in Table IIVI 





Triggers 


W 




TVa 




(n> 


2.567 GeV 
4.261 GeV 
5.754 GeV 


910 M 
1599 M 
5083 M 


1.6 - 2.15 GeV 
1.6 - 2.6 GeV 
1.6 ~ 2.6 GeV 


0.3 - 1.5 GeV2 
0.7 - 3.5 GeV^ 
1.3 - 5.4 GeV^ 


42000 
34000 
82000 


8000 
6500 
16000 


67% 
67% 
72% 



TABLE IV: Information regarding the different CLAS electroproduction data sets associated with this work, including the 
number of raw triggers, the W and extents of each data set, the number of hyperons (detected via the e'K'^p final state), 
and the average longitudinal polarization of the electron beam. Note that the 4.261 GeV data set sums together data acquired 
at beam energies of 4.056, 4.247, and 4.462 GeV. 

C. Particle identification 

The first level of event reconstruction required the identification of a viable electron candidate. This was done by 
requiring that a negatively charged particle - identified by its track curvature in the magnetic field of the spectrometer 
- be matched in time and space with hits in the SC, CC, and EC counters. A particle-tracking vertex cut was employed 
to ensure that the particle originated from the liquid- hydrogen target. In order to remove negatively charged pions 
from the electron candidate sample, a cut was placed on the ratio of the measured energy deposited in the fiducial 
region of the EC (accounting for the sampling fraction of the calorimeter) to the momentum of the particle. A further 
reduction in pion contamination was achieved by placing a minimum-ionizing cut on the energy measured in the EC. 

The first-level requirements for charged hadrons are that they have a track in the drift chamber and a matched, 
in-time hit in the SC in that same sector. For the final-state K'^ and p in this analysis, we require that the curvature 
for the K'^ and p tracks be consistent with a positively charged particle. We also require that the track originate 
from the target using a vertex cut. 

The algorithm used for hadron identification was slightly different between the elc and el-6 data sets. For the 
elc data set, the final-state particles were identified with momentum-dependent cuts on the momentum vs. mass 
distribution (to account for the worsening resolution of CLAS with increasing momentum). For the el-6 analysis, 
hadron identification was performed using a timing cut. The timing quantity of interest (dt ^ ti ~ t2) was the 
difference in the time between the measured flight time for a particle from the event vertex to the SC system (ii) 
and that expected for a given hadron type (^2)- The quantity ^2 was computed for all positively charged particles 
assuming the mass of the pion, kaon, and proton. A small 5t indicates that the correct mass hypothesis has been 



22 

made. The timing cuts are defined such that only one mass hypothesis can be satisfied for a given hadron. Fig. [5] 
shows the St plots used to identify final-state 7r+, A'+, and p candidates. 




fit 7T VS. P fit K VS. P fit p VS. P 

FIG. 5: Time diflFerence St (ns) between the SC system and the flight time calculated for specific hadron-mass hypotheses 
plotted against the hadron momentum P (GeV) for (a) tt'^ , (b) K'^ , and (c) p at 5.754 GeV. Each plot shows the timing cuts 
imposed upon the spectra for the given mass hypothesis. 

Prior to imposition of final particle identification cuts and in order to ensure an optimal resolution for the hyperon 
missing mass spectrum, the reconstructed momenta for the electrons and charged hadrons in the final state were 
corrected for small imperfections in the torus magnetic field map and the drift chamber alignment by using reactions 
with over-determined kinematics. The size of the momentum corrections {5p/p) for each of the final-state particles is 
on the order of 1%. 

D. Acceptance corrections 

In order to correct the yields for the detector acceptance, it is necessary to employ cuts that define the regions 
of CLAS where the detection efficiency is reasonably large and uniform. These fiducial cuts for both electrons and 
positive hadrons depend on momentum, angle, and torus field setting. For the electron, the CLAS acceptance is 
determined mostly by the limits of the azimuthal angle 4>e acceptance in each sector. The (/)e limits are determined 
by a marked drop in the collection efficiency of the CC at the edges of the detector. Additional fiducial cuts for all 
charged particles are designed to exclude regions of non-uniform acceptance from attenuation due to interactions with 
the mini-torus coils, the torus cryostat, or from the edges of the drift chamber acceptance. 

The acceptance correction was based on an analytic calculation that determined the geometrical acceptance factor 
on an event-by-event basis given the (j) acceptance of each final-state particle within the defined geometrical fiducial 
region. This factor accounted for losses due to kaon decays in-flight, bad scintillator paddles in the SC system, and 
the A pn^ branching ratio. Typical acceptances for the KY reactions requiring detection of the e'K^p flnal state 
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are at the level of 5% to 20%. A detailed comparison between the nominal geometric acceptance correction and a full 
GEANT acceptance function was performed for the elc analysis (see Ref. j^sl). Both methods were shown to have 
very similar functional forms. However, the beauty of the asymmetry approach employed for this analysis is that the 
results are relatively insensitive to the acceptance correction. Thus a much simpler analytic form was chosen over a 
full Monte Carlo approach. 



E. Hyperon yield extraction 

The reactions of interest are identified from missing-mass [MM) reconstructions of the e'A'+ final state. Shown 
in Fig. [SJa) is the MM{e'K^) distribution for the e'K^p final state at 5.754 GeV. This spectrum shows substantial, 
well-separated peaks for the ground state A and TP hyperons. The width of the A peak in this spectrum, summed 
over all and W , is about 11 MeV. Fig.lHI^b) shows the MNP{e' K^p) (missing mass squared) distribution. Here the 
final-state proton can come from the decay of the A(1115) (missing 7r~), the YP{\1%2) (missing 7r~7), or the A(1520) 
(missing K^). Fig. Elja) requires a cut on the MM'^{e' K^p) spectrum in the range from 0.007 to 0.065 GeV^, as 
shown in the correlation plot of Fig. [HIc), to reduce the contributions of particle misidentification background. The 
final K^K and A'+S° yields are then extracted through the fitting procedures described below. 




1.1 1.2 0.2 0.4 



MM (eK^)p MM' (eK>) MM' vs. MM 

FIG. 6: (Color online) e'K^p final state data from runs at 5.754 GeV. (a) Missing mass for p{e,e' K'^)X . (b) Missing mass 
squared for p{e,e'K+p)X. (c) MM^{e'K+p) vs. MM{e'K+) showing the cuts used to select the K+A and K+E° event 
samples. (Units in GeV.) 

The three components to the hyperon missing mass spectrum are the K^A final-state events, the ii'+E*' final-state 
events, and the particle- misidentification background (dominated by pions misidentified as kaons). As discussed in 
Section [iVi these individual contributions must be determined to extract the A and T,^ polarizations, and they have 
been measured through fits to the hyperon mass distributions. In this procedure, the A and peaks were fit using 
templates derived from a phase-space GEANT Monte Carlo simulation. The templates were generated with radiative 
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MM (eK*)p MM (eK*)p 

FIG. 7: (Color online) Sample hyperon spectrum fit results using hyperon templates derived from Monte Carlo (A:green curve, 
E'':red curve) and a third-order polynomial for the background (magenta curve). The smoothed templates have been allowed 
to shift along the mass axis and convoluted with a Gaussian to give the best per degree of freedom (x^) values for the fits. 
These distributions at 5.754 GeV are summed over all W and cos S^"*' for central values as indicated. The blue curve shows 
the full fit result. 



effects turned on, which is necessary to account for the A radiative tail beneath the peak. The background 
contributions in each bin were studied employing two different procedures. In the first, a background spectrum was 
derived from Monte Carlo using a phase space generator for multi-pion final states. The final-state 7r+ were then 
assigned the mass. The resultant MM(e'7r+) spectra were then sorted into the different analysis bins in Q^, W ^ 
and cos 6'^™-. The second approach employed a third-order polynomial to fit the backgrounds in the spectra. These 
two models for the background gave consistent answers, however the polynomial model was employed for the final fits 
as the Monte Carlo background distributions were statistically limited. 
The form of the hyperon spectrum fit in each analysis bin was given by 



MM = A ■ Ktemplate + B ■ 

^template + P(3)6efe, (31) 



where ^.template and "^template are the simulated hyperon distributions with weighting factors A and -B, respectively, 
and P(3)6cfe is a third-order polynomial describing the background. In performing these fits, the A and Yf Monte 
Carlo templates were allowed to shift up to ±10 MeV to match the data. In addition, the hyperon templates were 
individually convoluted with a Gaussian with a width chosen to minimize the of the fits in each bin. This was 
necessary as the resolution of the Monte Carlo was not a perfect match to the real data. Finally, the Monte Carlo 
templates were smoothed using a spline fit to remove the effects of statistical fluctuations in the simulation samples. 
Fits for two representative bins are shown in Fig. [71 
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FIG. 8: Ratio of the number of background counts to A counts in the A analysis mass window weighted by the uncertainty in 
the ratio for each of the kinematic bins in this analysis at 5.754 GeV. 
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FIG. 9: Ratios of the number of counts in the E'' analysis mass window for the data at 5.754 GeV. (a) TP / K ratio, (b) 
Background/S'' ratio. 



Fig.[8]shows the results of the yield fits in terms of F^ck — ^fccfc/^A for events in our A mass window {MM{e'K^) 
from 1.080 to 1.160 GeV) for each of our data sorts. The plot shows the distribution of the background ratio for each 
bin weighted by the uncertainty in the ratio. As the particle misidentification background was found to be relatively 
independent of kinematics, a single value of Fbcfc=3.3% has been employed for all analysis bins based on the weighted 
mean of Fig.[51 Additionally, from the spectrum fits we have extracted the ratios Fy = N^o /N\ and F^ck ~ Nhck/N^o 
in our mass window {MM{e'K^) from 1.175 to 1.213 GeV) (see Fig. [9]). These ratios are relatively independent 
of the kinematics. However, due to the sensitivity of to the number of A events in the E'^ mass window, we have 
employed the measured E*^, A, and background yields in the analysis for each kinematic bin. 
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F. Background polarization corrections 

Once the number of A, T,'^, and background events are determined in each of the respective hyperon mass windows, 
the measured polarization must be corrected as discussed in Sec. IIV CI The background polarization was measured 
by sorting data from the ep — > e'n^pX reaction and assigning the mass to the tt^ events. The analysis procedure 
then followed all of the same steps and procedures as for the hyperon polarization analysis. The measured background 
polarization for one typical data sort is shown in Fig.[Tn]as a function ofW (summed over all other kinematic variables). 
The results are consistent with Vl^j.—Q for all axes in all of the sorts investigated. 

Within the tight A cuts, there was no measurable level of contamination. However, within the mass window 
(1.175 to 1.213 GeV) there is significant contamination from both the A radiative tail and pion background (see 
Table \V\ . To correct the E*' polarization for the A tail, Vj^ was determined following our nominal prescription for 
determining the transferred A polarization, where the A data were binned in the same bins as the S° data. This value 
of V'jy was then used in Eq. (P5)) to calculate the corrected S° polarization. 



W Bin (GeV) 




Na 




1.825 


3067±79 


2746±20 


92±27 


1.975 


3872±85 


2544±22 


214±28 


2.125 


2425±72 


2074±20 


329±24 


2.275 


3160±80 


2050±20 


501±26 


2.470 


3413±80 


1758±19 


408±19 



TABLE V: Results from the 5.754 GeV yield fits for TVs, A'^a, and Nbck with statistical uncertainties in a ±2a E mass window 
for a sort dividing the S" analysis into 5 bins in W. 

G. Radiative corrections 

No radiative corrections have been applied to the data in this analysis. These have purposefully been avoided 
by employing relatively tight cuts on the reconstructed hyperon spectrum for the K'^A and K'^T,'^ events and by 
accounting for the K^A radiative tail events within the if+E^ event sample. This is expected to be a reasonable 
approach as the radiative effects are independent of the beam helicity and thus should effectively cancel out of the 
asymmetry calculation. With our relatively tight hyperon mass cuts, the maximum radiated photon energy is only 
about 50 MeV, which has a negligible impact on our computed cos 9?^^ values with respect to each quantization axis. 
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FIG. 10: Polarization V' vs. W for the pion misidentification background in the {x',y',z') coordinate system for the data 
summed over aU Q2 and dn^"'. 

H. Bin averaging corrections 

The bin sizes for W, Q^, and cos0|^>™' in this analysis were chosen to roughly equalize the statistical precision 
of each polarization data point. To account for the finite bin sizes and the variation of the cross section over the 
bins, we quote our polarization results at the bin means. The bin mean was determined by measuring the mean of 
the acceptance-corrected yield distribution over the kinematic bins of interest. As might be expected, the largest 
differences between the bin mean and the bin center occur where the bins are larger. The kinematic bin means for 
each data sort are given in Section IVlII 



VI. SYSTEMATIC UNCERTAINTY ANALYSIS 

In this section we examine the sources of systematic uncertainty that affect the extracted polarization observables 
for the 5.754 GeV data set. The assigned systematics for the 4.261 GeV K^A data set are described in Ref. 49 1 



and are given by SV'sys < 0.084. The contributions to the total systematic uncertainty belong to one of four general 
categories: Polarization extraction, beam-related factors, acceptance function, and background contributions. As 
the statistics for the K^A final state dominate those for the i^T+E" final state, the determination of the assigned 
systematics for both final states is based on analysis of the if+A data. The exception to this is the assignment of a 
separate systematic for the background contributions. The final systematic uncertainty compilation for the 5.754 GeV 
measurements is given in Table IVII 
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The procedure used to assign a systematic uncertainty to each source within a given category is to compare 
the measured polarization V' for all kinematic bins with the nominal analysis cuts or procedures (nom) to that 
with modified cuts or procedures (mod). The spread in the difference of the polarization over all data points, 
AT" — 'P'lom ~ 'P'mod^ used as a measure of the systematic uncertainty for a given source. The estimated common 
uncertainty for a given source is the weighted root-mean-square (r.m.s.) of AV' for all points given by 



' sys 



where the sums are over all N data points and 5Vl is the statistical uncertainty of the i*'' data point. In the studies 
done for this analysis, a common systematic uncertainty is applied for all data points as the kinematic dependence 
of the AV' distributions was found to be minimal. In each of the systematic uncertainty studies performed for this 
analysis, the widths of the AV' distributions were much larger than the measured centroids, which are all consistent 
with zero. Thus the assignments are believed to be rather conservative. 

A. Polarization extraction 

The polarization has been determined by using two different analysis approaches. The nominal technique is the 
asymmetry approach described in Section IIVI An alternative approach is to extract the polarization from the ratio 
of the acceptance-corrected, helicity-gated yields via 

The difference between these two techniques resulted in an estimated systematic uncertainty of SV'syg = 0.008. 

A systematic uncertainty arises from the somewhat arbitrary choice made for the cos bin size. Nominally the 
data were sorted into six bins in the rest frame proton angle. A comparison of the nominal polarization results with 
the extraction from a sort with four and eight bins in this variable resulted in a weighted r.m.s. of SV'sys = 0.018. 
The difference in the polarization results is effectively due to the fitting algorithm employed in which the centroids of 
the cos^^^ bins are assigned to the center of the bin. When the number of bins is reduced, the fit results are more 
sensitive to the bin content. 

The final systematic uncertainty contribution in this category arises due to the uncertainty in the weak decay 
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asymmetry parameter Q!a. This uncertainty gives rise to a scale- type uncertainty on the extracted polarization (the 
same for both A and T,^ hyperons) given by SV^y^ = \VY\Sa\/aA = 0.02\Vy\- 

B. Beam-related factors 

There are two possible contributions to the systematic uncertainty related to the beam. The first factor is associated 
with the beam polarization measurement from the M0ller polarimeter system. This arises from the uncertainty in 
the M0ller target foil polarization, the statistical uncertainty in the measurements, as well as a contribution from 
variations of the polarization measurements over time. These contributions have been estimated for CLAS polarization 
measurements to be 4%. The associated uncertainty in the hyperon polarization is SP^y^ = IPyl^Pb/Pb = 0.04|'Py|. 

The second beam-related contribution is the beam charge asymmetry that results from a difference in the electron 
beam intensity for the two beam helicity states. Prom studies of the 5.754 GeV data set, the beam charge asymmetry 
was below the 10"'^ level and no detectable difference between the helicity-gated live times was found, thus no 
systematic contribution was assigned. 

C. Acceptance function 

There are several factors that go into the systematic uncertainty associated with the form of our acceptance 
correction and with the choices made to implement this correction, which include the specific form of the fiducial 
cuts used to define the azimuthal extent of the acceptance as a function of polar angle and the minimum acceptance 
cutoff. In order to assign a systematic uncertainty associated with the acceptance correction, we have compared the 
extracted polarizations with and without the geometric acceptance corrections. The r.m.s. width of the difference 
distribution was assigned as the systematic uncertainty for the acceptance correction. This value, SP'^y^ = 0.033, is 
believed to be a very conservative estimate. 

To study the effects of the fiducial cuts employed to define the azimuthal acceptance for electrons and hadrons, two 
different sets of fiducial cuts were defined in the analysis. A loose cut (the nominal cut) was designed to define the 
azimuthal acceptance edge of CLAS as a function of momentum, and a second cut was designed to be several degrees 
tighter than the nominal cut. Comparisons of the extracted polarizations between these two cut definitions gave an 
r.m.s. width of SP'^yg = 0.020, which represents the assigned systematic uncertainty. 

The minimum acceptance cutoff translates into a maximum acceptance weight. The minimum acceptance cutoff 



30 

was nominally set at 10% (a somewhat arbitrary choice) for the e'K~^p final state. For our study, we varied the 
acceptance cutoff by ±20% relative to the nominal cutoff value. The assigned systematic uncertainty, given by the 
r.m.s. width of the polarization difference distribution, is SV^yg = 0.025. 

D. Background contributions 

1. e'K'^A final state 

Our analysis of the backgrounds found no measurable level of T,^ contamination within our final K~^A event sample. 
However, there is a few percent contamination of 7r+ misidentification events that remain beneath the A peak that 
serve to dilute the measured A polarization. To estimate the systematic uncertainty associated with our subtraction 
technique, we have compared our nominal polarization results to results obtained assuming no pion background. 
Clearly, this would result in an overestimate of the systematic uncertainty so we have used one-half of the difference, 
or SV'syg — 0.009. While an arbitrary choice, this represents a conservative estimate and is small compared to other 
sources of systematic uncertainties. 

S. e'K+YP final state 

The uncertainties in the backgrounds from K'^ K and pion misidentification have a much bigger impact on the 
extracted YP polarization compared to the A analysis. Therefore, it is important to study these effects separately 
for this final state. Our approach to assign a systematic uncertainty due to the fit uncertainties of the contributing 
backgrounds beneath the YP is to allow the extracted yields to vary by ±10% from the fit value and to study the 
effect on the extracted polarization. Variations of the background levels of ±10% amounted to variations on the 
fitted background yields of A'a ± 2cr7VA ^-nd Ni,ck i 2(7^^^^. Our studies indicated that the maximum change in the 
measured YP polarization was ±0.10, which we have assigned as the associated systematic uncertainty SV'gy^. 

E. Final systematic uncertainty accounting 

Our final systematic uncertainty accounting for the 5.754 GeV A and V' data is included in Table fVll listing 
all of the sources discussed above. The final value for the total systematic uncertainty results from adding all the 
individual contributions in quadrature. (Additions in quadrature in Table IVTl are represented by the notation ©). 
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Category 


Contribution 


Systematic Uncertainty 


Polcirization Extraction 


Functional Form 


U.UUo 




Bin Size 


U.Ulo 




Asymmetry Parameter 


U.UZ I-'y 


Beam-Related Factors 


Beam Polarization 


0.04 V!y 


Acceptance Function 


Fiducial Cut Form 


U.UzU 




Acceptance Correction 


U.UoO 




Acceptance Cutoff 


U.UZO 


Background Contributions 


Pion and A 
contamination 


0.009 (A), 0.100 {YP) 


( Total Systematic Uncertainty ) 


0.051 (A), 0.112 {TP) 






® 0.0457'(. 



TABLE VI: Summary table of the systematic uncertainty assignments SV'sya for the measured A and YP polarizations at 
5.754 GeV. Note that the contributions from the uncertainties in the beam polarization and weak decay asymmetry parameter 
(given by 0.045'Py) are added in quadrature (represented by the © notation) to the other sources. 



One way to verify the veracity of the final systematic uncertainty assignment is to look at the deviations of the 
normal components of the extracted A and YP polarizations (i.e. along the y' and y axes). Averaged over all analysis 
bins, the weighted mean of the V'yi and Vy components for the A is 0.067 and for the TP is 0.134. Both of these values 



are consistent with our total systematic uncertainty assignments in Table rVTl The extracted normal components for 
one of our data sorts for the A and TP hyperons are shown in Fig. 1111 

W=1.99GeV cos (Jk"" =0.43 W=2.09 GeV cos 0k "" = 0.36 




(GeV^) 



(GeV^) 



FIG. 11: (Color online) Distributions of the transferred A (a) and YP (b) polarization components V'yi vs. for two of our 
analysis bins. The error bars shown represent the statistical uncertainties only. 
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VII. RESULTS AND DISCUSSION 



A. A polarization transfer 



Our results for the transferred A polarization acquired at a beam energy of 5.754 GeV are shown in Figs . fT2l through 
1151 compared to several model calculations. The error bars in these figures include statistical but not systematic 
uncertainties, which we estimate to be 0.051 ® 0.0457^^ on the polarization. The full data set is contained in the 



CLAS database 



50|. 



Figs. [T^ and [T51 show the dependence of "P^, ^, and V'^ ^ with respect to cosdjf"- for the three bin-averaged W and 
values indicated in the figure. Fig. [T^] shows that the value of V'^, decreases smoothly with increasing scattering 
angle, whereas V'^r decreases with increasing angle until cos0|^'"' w 0.8, at which point it levels off to a value of about 
—0.5 over the range of cos6'^™' covered by the experiment. The fact that the V'^, data approach zero at cosf?^™' = 1 
is simply a result of angular- momentum conservation, which also requires V'^i = at cos0^™' = — 1. The dependence 
of the polarization along the {x, z) axes in Fig. [12] is qualitatively different. The polarization along V'^ is roughly 
zero everywhere, whereas V'^ is relatively constant (at least over the angle range cos 6'^™' > where the statistics are 
reasonable) with an average value of ~0.6. This may be hinting at a simple reaction mechanism (see Section fVlI Pp . 

The polarization with respect to W is shown in Fig.[T3]for a beam energy of 5.754 GeV. Note that there is a strong 
dependence of the bin-averaged cos0|^™- value with respect to W in these kinematics. The central cos^f^™- values 
extracted from the analysis are reasonably represented by a fit to a second-order polynomial in W , with the fit shown 
in Fig. [Ml In the plot of V'^, (Fig. [T4)) we see that the polarization rises steadily from zero near threshold, followed 
by a dip at around 1.9 — 2.0 GeV, and then remains constant at about 0.5 over the rest of the range. The V'^, data 
are relatively constant at about —0.5 over most of the W range. We also note that the magnitude of V'^, is nearly 
equal to V'^,, although with opposite sign, indicating equal strength in the i?f.j., and responses. With respect 

to the [x, z) axes (see Fig. I14p . the polarization is roughly 0.6 and relatively constant along z and is consistent with 
zero along x. This latter point indicates either a perfect cancellation of the Rlt' response functions (see Table [llT]) 
or that they are each nearly zero. 

The polarization as a function of is shown in Fig. [THl The data are rather featureless and indicate almost no 
dependence. Both V'^j and V'^ are roughly 0.5, while V'^, ~ and V'^ is consistent with zero. 

The V'l^ data in Figs. [T^ through [T51 are compared against the theoretical models introduced in Section [Til The 
MB hadrodynamic model 42] is indicated by the solid-black lines, the GLV Regge model 51l| is indicated by the 
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FIG. 12: (Color online) Translerred A polarization components V' with respect to the (a;', 2') axes vs. cos^J 
averaged W jQ^ values as indicated for a beam energy of 5.754 GeV. The curves are calculations from the MB isobar model 
[131 (solid - black), the GLV Regge model [IJ] (short dash - blue), and the RPR model [s^l variant including a Pi 1 (1900) state 
(dot-dash - red) and a Z)i3(1900) state (long dash - black). 




FIG. 13: (Color online) Transferred A polarization components V' with respect to the (a;, z) axes vs. cos S^™' for three bin- 
averaged W IQ^ values as indicated for a beam energy of 5.754 GeV. The model calculations are as indicated in Fig. 1121 
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FIG. 14: (Color online) Translerred A polarization components V' with respect to the {x' , z') axes (upper left panels) and {x, z) 
axes (upper right panels) vs. W (GeV) for ((3^)=2.54 GeV^ for a beam energy of 5.754 GeV. The bin-averaged cosS^™' values 
for each W point are shown in the lower plot. A fit to {cos 9%"^') is provided by a second-order polynomial in W as indicated 
on the plot. The model calculations are as indicated in Fig. 1121 



short dash-blue lines, and the RPR model 52] is indicated by the dot-dash-red lines (Pn model variant) and by 
the long-dash-black lines {Di^ model variant) . The calculations qualitatively match the sign and trends of the data, 
but detailed comparisons indicate that these new polarization data can be used to further tune the models (e.g. the 
resonance parameters in the MB model and the RPR model) or indicate shortcomings in the dynamical description 
of the data (i.e. the pure t-channel description of the GLV Regge model). 

Detailed comparisons of the individual models to these data are also useful to indicate specific shortcomings of 
the models. For example, comparisons of the MB model to the data show problems with the parameters for the 
resonances included below 2.0 GeV as indicated by both V'^, and V'^,. The models also mostly fail to reproduce the 
data as a function of (Fig. fT5| . which could indicate problems with the modeling of the non- resonant strength 
with increasing or the description of the evolution of the hadronic form factors. 
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(GeV^) 
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FIG. 15: (Color online) Transferred A polarization components V' with respect to the {x ,z') axes (left panels) and {x,z) axes 
(right panels) vs. (GeV^) for (W) = 1.99 GeV and (cos Sx™ ) =0.43 for a beam energy of 5.754 GeV. The model calculations 
are as indicated in Fig. 1121 



Comparisons of the GLV Regge model to the data indicate that a purely i-channel description of the K reaction 
is not adequate to reproduce the polarization results, even for this channel suspected to be predominantly governed by 

n 

K and K* exchange ,6]. The s-channel resonance contributions still have important consequences for the interference 
observables. Whereas the GLV model produces a very smooth behavior for V' vs. W , cos6'|^™', and , the model 
typically underpredicts the strength and does not account for the detailed trends in the data. In some cases (e.g. 
Figs. [Tallin and [T5|) . the GLV model has the wrong sign compared to the data or has the wrong slope. 

For the RPR calculation, two model variants are compared with the data. One employs a Pii(1900) state. As was 



also seen in comparison with the CLAS (Jlt' data [Sj], the model variant with the Pn (dot-dash - red) is strongly 
ruled out by the polarization data. The second model variant (long dash - black) employs the Z?i3(1900) state 
proposed by Mart and Bennhold (see Section [Tl| . This model provides a reasonable description of the polarization 
data over the full kinematic phase space. The only issue with this model, which cannot be fully clarified by these 
data, is the strong interference effects seen in the calculations at higher W (see Figs. [T^ and [T5)) . 



B. YP polarization transfer 

Our results for the transferred S*^ polarization acquired at a beam energy of 5.754 GeV are shown in Figs. [TC] through 
[TSl compared to several model calculations. The error bars in these figures include statistical but not systematic 
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uncertainties, which we estimate to be 0.112© 0.045'P^ on the polarization (see Section |Vl|. The full data set is 
contained in the CLAS database 5C | . 

Fig. [TBI shows the dependence of V'^, ^, and V'^ ^ with respect to cos6'^™' and Fig. [T7l shows the polarization with 
respect to W . As with the A results, cos 6'^™' values at each W point are well represented by a low-order polynomial 
in W (given in Fig. [T7|. Fig. [TSl shows the polarization with respect to . 
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FIG. 16: (Color online) Transferred S° polarization components V' with respect to the {x',z') (left) and {x,z) (right) axes 
vs. cosO'jf"' at (Q^)=2.5 GeV'^ and {W)—2.1 GeV for a beam energy of 5.754 GeV. The curves are calculations from the MB 
isobar model [42| (solid - black), the GLV Regge model 51] (short dash - blue), and the RPR model [52i] with a missing 
Pii(1900) state (dot-dash - red) and a missing Z)i3(1900) state (long dash - black). 



An important point to note is the rather sizable statistical uncertainties on the Yp data. This arises not due to the 
limitations of the S]° data sample (which has roughly 3000 counts in each analysis bin), but rather due to the scaling 
by QfA/as in Ea. (|26p . It should also be made clear that the el-6 running period at CLAS represented a very lengthy 
run period (extending over 4 months) that recorded nearly 5 billion triggers. Given the effect of the A background 
beneath the Yl^ peak, CLAS will likely not be able to provide more precise electroproduction data for this observable. 

Figs. [TBI through [TSl indicate that the measured transferred polarization tends to have the same sign as the A 
polarization for the different axes. More quantitative statements cannot be made given the statistical quality of the 
data. In comparing the data to the models, with the possible exception of the MB model (solid black), one sees that 
they qualitatively match the sign and trends of the data, as was the case for the A data. However, as the V'y^ observable 
has not been measured before, the results can still serve to provide at least loose constraints on the theoretical models. 
Certainly the V'j^ data can be used to improve the knowledge of the contributing N* states by allowing for improved 
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FIG. 17: (Color online) Transferred E° polarization components P' with respect to the {x' , z') (left) and {x, z) (right) axes vs. W 
(GeV) at (Q^)=2.5 GeV^ for a beam energy of 5.754 GeV. The bin-averaged cos ^^^^ value for each W point is represented 
by a second-order polynomial in W . The model calculations are as indicated in Fig. 1161 



descriptions of the associated form factors and the qkna coupling strengths. These improvements can then be used 
in the S° modeling, which typically employ the same set of N* states for both the K^A and X+S*^ final states. 



C. Extraction of the ratio of (Tl/tt 



In order to extract Ra = (JlI^t at cos6'^'"' = 1, we must first extrapolate V'^j (or V'^ to cos6'|^'"' = 1. However, 
because of statistical fiuctuations in the data and finite angle resolution effects, extrapolations for V'^, give slightly 



different results than extrapolations for V'^ . Following the procedure defined in Ref. 



23| , the extrapolation is actually 



performed by summing the V'^, and V'^ components into a new quantity Rsum (see Table Hill for component definitions) 
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FIG. 18: (Color online) Transferred E'' polarization components V' with respect to the (x! ,z'^ (left) and (a;, z) (right) axes 
vs. (GeV^) at (W)=2.09 GeV and (cosS^™ ) — 0.36 at a beam energy at 5.754 GeV. The model calculations are as indicated 
in Fig. [H 



given by 



-Rsum — ^ ^ = /C[(l 

Co 



)R 



z'O 
TT' 



jDx'O ■ nc.m.-\ 



(34) 



From the extrapolated value of Rsum at cos0J>™- = 1, Eq. (IM)l can then be inverted to determine V'^, (or equivalently 
V'^), which in turn is used to extract using Eq.Q. 

An additional benefit of using this form is that Eqs.(I7]) and ^ provide important and useful constraints on Rsum 
aX X = cosOj^™- ~ ±1. At X ~ —1, the sum of the polarizations must be zero, according to Eq.(l7]), leading to 
Rsumix — —1) — 0. Since both ul and ctt must be positive definite, then R^ (see Eq.®) must also be positive 
definite, which leads to Rsum < 2(Ty for x = 1. Note also that as both o-^ and ax must be positive definite, Eq.® 
constrains both the A and E° hyperon polarization components V'^ to be between and cq. This argument confirms 
the sign of the hyperon polarization is correct in our analysis. 

Besides the explicit 0^™- dependence shown in Ea. (p4)) and in the response functions, the CGLN amplitudes contain 



35| . This suggests that Eq. ([M|) can then be fit with 



additional 6'^™- dependence (as well as and W dependence) 
polynomials in a:: = cos6'^™', provided we have prior knowledge of the au term. In the case of the 4.261 GeV data, 
we can use the previously published CLAS results 6], while for the 5.754 GeV data, we have to use models to provide 
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ajj at our kinematic points. 

The number of terms to include in a polynomial fit to Eq. (|34p is ultimately governed by the reaction dynamics. 
The explicit dependence alone suggests at least a third-order polynomial. However, given the limited number of 
polarization data points, the number of terms in any fit leading to a meaningful extrapolation to cos0|^™' = 1 must 
also be limited. We begin by considering third-order fits of the form 

Rsum = ao + aix + a^x^ + a:iX^ , (35) 

where 0^=0^3 represent the fit coefficients. However, applying the constraint Rsum = at x = — 1 implies ai = 
ao + 02 - 03. 

We have done a series of fits to the data points representing Rsum in which we varied the number of terms in the 
fits, while imposing a penalty on the if a fit returned an unphysical value at a; = 1. The penalty was chosen to 
be large enough to force non-negative values of R^- It should be noted that only one of the fits (5.754 GeV data 
at W=1.75 GeV and Q^=2.61 GeV^) required the imposition of a penalty. In determining the optimal number of 
parameters in the fit for each W, we simply used the number of parameters that produced the smallest minimized xt 
(x^ per degree of freedom). All three of the 4.261 GeV fits favored a second-order fit (03 = 0), while all three of the 
5.754 GeV fits favored a third-order fit. 

The V'^, and V'^ data from the 4.261 GeV data set are shown in Fig. [12] with respect to cos0|^™'. The full set of 
4.261 GeV polarization transfer data for the {x',z') and {x,z) axes with respect to W and cos^^™- is provided in 
Ref. 50]. The V'^, and V'^ data from the 5.754 GeV data set are shown in Figs. [T^] and [T51 We should point out 
that the V'^, and V'^ results come from the same data. Therefore, these observables are not independent. They do, 
however, measure different quantities (as seen in Table IhB since they are projections onto different axes. In adding 
these together to form Rsumi the uncertainties from V'^, and V'^ were added together. 

The results of our fits to the 4.261 GeV and 5.754 GeV data are shown in Fig. [^D] (heavy solid lines) along with 
an error band (light solid lines). The error bands include uncertainties both from the fitting of Eq. (|35p . and, for 
the 4.261 GeV data, contributions from uncertainties in the fits of the cross section data. The latter contribution to 
the uncertainties is about half that of the former. The error band indicates that the extrapolation to x = 1 is well 
constrained. For the 5.754 GeV data we display the fit using the calculated cross section for one particular choice of 
the MB model 26], which allows for different choices of form factors and couplings. 

Table IVTIl shows the resulting x^, the polarization extrapolated to a; = 1, and R^- Since the 5.754 GeV data 
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FIG. 19: (Color online) Transferred A polarization components V' with respect to the {z',z) axes vs. cosS^™' for three 
bin-averaged W/Q^ values as indicated for a beam energy of 4.261 GeV. 



Eb (GeV) 


(W) GeV 


(g2) GeV2 


xl 




Ra 




1.72 


1.63 


0.77 


0.451±0.066 


0.533±0.270±0.326 


4.261 


1.89 


1.58 


5.69 


0.440±0.063 


0.870±0.329±0.401 




2.18 


1.45 


0.87 


0.486±0.062 


1.348±0.404±0.515 




1.75 


2.61 


1.11 


0.607±0.070 


0.000±0.092±0.156 


5.754 


1.98 


2.56 


2.76 


0.610±0.065 


0.176±0.088±0.209 




2.31 


2.31 


2.79 


0.470±0.053 


0.637±0.120±0.445 



TABLE VII: Transferred polarization at a; = cosO'jf"' — 1.0 extrapolated from the fits described in the text along with 
the resulting value of the ratio of longitudinal to transverse structure functions. Uncertainties on V'^, ^ are the combined 
uncertainties arising from the fit to the polarization data and the uncertainties in the cross section data. The first uncertainty 
on Ra is the statistical uncertainty (from the fit), while the second represents an estimated systematic uncertainty. 



required a model for au, we repeated the fit for five different parameter set choices within the framework of the MB 
model Thus, the 5.754 GeV results in the table reflect the average values of x^, the polarization extrapolated to 
X — I, and for different models. We estimated the model uncertainty by using the standard deviation of Ra from 
using the five different models. 

Inserting the extrapolated polarizations into Eq.([9]), we can determine the ratio Ra-. These values are shown in 
the last column of Table IVIIl along with the combined uncertainties of the polarization and cross section fits, and 
an estimated systematic uncertainty. The systematic uncertainty includes a contribution assuming a 10% relative 
systematic uncertainty in the polarization data. 
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4.261 GeV Data 



5.754 GeV Data 




FIG. 20: (Color online) Rsu^ (defined in Eq.JHll) vs. cos6'ff™- for the 4.261 GeV and 5.754 GeV data along with our fits 
(heavy solid lines) and the error band resulting from the fit uncertainties (light solid lines). The dashed red line in the upper 
right panel indicates the result of removing the x = 1 constraint in the fit. 



The resulting values for R„ are plotted in Fig. [^T] For comparison, we have also included the previously published 



data a, |2J, |37[. However, only the filled points are at or near the common value of W «1.84 GeV. Other than these 



filled points, one should not take any trends in the data too seriously since the data from this analysis cover a large 
range in W (1.72 to 2.31 GeV). 
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FIG. 21: (Color online) Ratio of longitudinal to transverse structure functions at — 0° vs. Q^. The blue cross data points 

(about Q^=1.5 GeV'^) are for 4.261 GeV and the red cross points (about Q^=2.5 GeV^) are for 5.754 GeV. The inner error 
bars on our points represent the statistical uncertainties arising from the fit and the outer error bars represent the combination 
of statistical and estimated systematic uncertainties. The triangles are data from CLAS !6|, the solid circles are data from 
Mohring et al. [53], the solid squares are data from Bebek et al. [37j . and the diamonds are data from Raue and Carman [23l |. 
All of the filled points are near a common value of W «1.84 GeV. 



Our new results from the data sets at 4.261 and 5.754 GeV are in reasonable accord with the existing measurements 
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of (Tl /ct from Bebek et al. [37'| and Mohring et al. 



24j . Looking at the results for W ftil.84 GeV, the ratio rises with 
up to ^ 1.5 GeV^ and thereafter seems to fah off, suggesting an interesting and non-trivial dependence on Q^, 
but the measurement accuracy is not adequate to quantify this observation. Alternatively, we point out that our data 
also seem to suggest a rapid rise of ctl/ct with W as was suggested in our previous publication 6|, however again, 
we lack the statistical and systematic precison to make a more definitive conclusion. Note that the aL/cT data from 
Ref. [6| at Q^ — l.Q GeV^ cannot be directly compared to these data as the most forward angle point in that work is 
cos6iff™=0.90. 

The data of Fig. [21] imply that for at least a limited interval, the longitudinal structure function becomes sizable. 



This structure function is expected to be very sensitive to the kaon form factor 



53| . A recently conducted experiment 



in Hall A at Jefferson Laboratory 54] has as one of its main goals a Rosenbluth separation at several values of 



55j of the kaon form factor. However, this method relies 



momentum transfer t leading to a Chew-Low extrapolation 
on having small relative uncertainties for ct^, which will not be the case when is itself small. These new results 
indicate that the successful extraction of the form factor may only be possible in a limited kinematic range. 



D. Partonic models of the process 

All of the models introduced thus far in this work have been used to indicate the strong sensitivity of these 
polarization data to the underlying s-channel resonant terms that contribute in the intermediate state of the — > 
K^Y process. The precision and broad kinematic coverage of the data from Ref. [6^ have indicated that the K^A final 
state is dominated by ^-channel kaon exchange. However, there are important contributions from s-channel processes 
that must be taken into account to describe both the cross section and polarization data in detail. 

In contradistinction to the hadronic models, and as noted earlier and introduced in Ref. 10|, our data indicate the 



A polarization is maximal along the virtual photon direction (see results for in Figs. [T51 1141 and llSp. suggesting 
a simple phenomenology. In fact, the A polarization is essentially unity if the virtual photon depolarization factor 
is taken into account (see Section flVDp . The lack of a strong W and dependence is an indication that the data 
might be more economically described in a flux-tube strong-decay framework. There is growing evidence that the 
relevant degrees of freedom to describe the phenomenology of hadronic decays are constituent quarks bound by a 
gluonic flux-tube 



56[. Properties of the flux-tube can be determined by studying qq pair production, since this is 
widely believed to produce the color field neutralization that breaks the flux-tube. Since the 1970's, it has been argued 
that a quark pair with vacuum quantum numbers is responsible for breaking the color flux-tube (the ^Pq model ,57j). 



43 






d 






u 


U 


u 





b) 




FIG. 22: (a) A model of the reaction where a circularly polarized virtual photon strikes an oppositely polarized u quark inside 
the proton. The spin of this quark flips and the quark recoils from its neighbors. An ss quark pair is created from a, = 0~ 
two-gluon exchange (in lowest order) to produce the final-state and A hyperon. (b) A model of the reaction where an ss 
quark pair is produced from a circularly polarized real photon that hadronizes such that the s quark in the A retains its full 
polarization after being "precessed" by a spin-orbit interaction, while the s quark ends up in the spinless kaon. In both pictures 
the shaded band represents a spinless u ~ d diquark system. 



This simple phenomenology of the A polarization data has led two groups within the CLAS Collaboration to 
develop semi-classical models based on partonic degrees of freedom to describe the associated reaction mechanism. 
In the model of Carman et al. 
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58| (shown in Fig. [221 all, it is assumed that the cross section is dominated by 
photoabsorption by a m quark. Due to the helicity-conserving vector interaction, the u quark becomes polarized along 
the photon direction {+z). Hadronization into the K^A final state proceeds with the production of an ss pair that 
breaks the color flux-tube. Because the u quark hadronizes as a pseudoscalar , the s quark spin is required to 
be opposite to that of the u quark, i.e. in the —z direction. In the non-relativistic quark model, the entire spin of 
the A is carried by the s quark. Since the A polarization is in the +z direction, seen by the fact that V'^ > 0, it 
was concluded that the s and s spins were anti-aligned when they were created, if the hadronization process did not 
flip or rotate their spins. Note that Liang and Boros also posit a two-step process for the production of transversely 
polarized A hyperons in the exclusive pp pK^A reaction 59|, and come to a similar conclusion that the ss quark 



pair must have been produced with spins anti-aligned. More recently, A polarization has been interpreted within an 
instanton interaction model jsol, which also is assumed to occur with the production of an anti-aligned ss quark pair 



25|, 



A dominance of spin ant i- alignment for the s and s quarks would not be consistent with the S = 1 '^Pq operator 
which predicts a 2:1 mixture of ss quarks produced with spins aligned vs. anti-aligned if the orbital substates are 
equally populated. Along with other observations of failure of the ^Pq model (e.g. explaining 7r2 ptu decay |25|), 
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the applicability of the ^Po model in describing all hadronic decays is brought into doubt if this model is appropriate. 

Extensive photoproduction data for the transferred polarization for the K final state has also been published 
from CLAS 9|. These data also indicate that the A polarization is predominantly in the direction of the spin of 
the incoming photon, independent of the center-of-mass energy or meson production angle. Based on these data, 
Schumacher has introduced a different model 
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6l| to explain the A polarization results. In this model, shown in 
Fig-llUb), the produced ss pair is created in a '^S'l configuration (J=l, 5=1, L=0, i.e. J'^ = 1^). Here, following the 
principle of vector meson dominance, the real photon fluctuates into a virtual <j) meson that carries the polarization of 
the incident photon. Therefore, the quark spins are in the direction of the spin of the photon before the hadronization 
interaction. The s quark of the pair merges with the unpolarized diquark within the target proton to form the A 
baryon, and the s quark merges with the remnant u quark of the proton to form a spinless meson. 

The two model interpretations, while able to predict the correct sign for the A polarization transfer, nevertheless 
describe very different physical processes. Both assume that the mechanism of spin transfer to the A hyperon involves 
a spectator = 0+ diquark system. The main difference is the role of the third quark. Neither model specifies a 
detailed dynamical mechanism. If we take the gluonic degrees of freedom into consideration, the model of Carman et 



al. lOj can be realized in terms of a possible mechanism in which a colorless J'" — 0^ two-gluon subsystem is emitted 
from the spectator diquark system and produces the ss pair as illustrated in Fig. I22r a) . To the same order of gluon 
coupling, the model of Schumacher 



6lj is the quark-exchange mechanism illustrated in Fig. I22r b) , again mediated by 
a two-gluon exchange. The amplitudes corresponding to these diagrams may both be present in the production, in 
principle, and could contribute at different levels depending on the reaction kinematics. 

Extending these studies to the X*+A exclusive final state should be revealing. In the Carman et al. model, the spin 
of the u quark is unchanged when changing from a scalar K'^ to a vector K*^ . If the ss quark pair is produced with 
spins anti-aligned, then the spin direction of the A should flip. On the other hand, in the Schumacher model, the u 
quark in the kaon is only a spectator; changing its spin direction - and thus changing the to a K*^ - should not 
change the A spin direction. Thus there are ways to disentangle the relative contributions and to better understand 
the reaction mechanism and dynamics underlying the associated strangeness production reaction. Analyses at CLAS 
are underway to extract the polarization transfer to the hyperon in the A'*+A final state. 

In developing the quark model interpretations of polarization transfer, we also need to consider the phenomenology 
of the K^YP results. As shown in this work, and much more clearly in the CLAS photoproduction data 9], the AT+E" 
polarization transfer is very similar in magnitude and sign to the A'+A data. We might expect that when the A and 
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polarization transfers in these reactions point in opposite directions and have the same magnitudes, this would 
then give more weight to the modeling of the polarizations originating from a quark level interaction, in particular 
associated with the strange quark spin. However, even though this is not what is observed in the data, we should 
realize that the spin state of the hyperon is not determined by the strange quark alone, but a combination of the s 
quark spin and the triplet ud quark spin. Thus the models of Fig.tSHa-re not directly applicable for X+S^ production. 

Understanding a process of this sort through partonic models can shed light on quark-gluon dynamics in a domain 
usually thought to be dominated by traditional meson and baryon degrees of freedom. These issues are relevant to 
better understand strong interactions and hadroproduction in general due to the non-perturbative nature of QCD 
for CLAS kinematics. We eagerly await further experimental studies and new theoretical efforts to understand which 
multi-giuonic degrees of freedom dominate in quark-pair creation and their role in strangeness production, as well as 
the appropriate mechanism (or mechanisms) for the dynamics of spin transfer in hyperon production. 



VIII. SUMMARY AND CONCLUSIONS 

In this paper we have provided extensive new data at 4.261 and 5.754 GeV for the beam-recoil hyperon polarization 
transfer for the reaction p{e, e'K^)A studying its dependence on the kinematic variables Q^, W, and cos0^™'. These 



data add to the earlier 2.567 GeV CLAS data results from Ref. In addition, we have provided the first-ever 

polarization transfer data for the reaction p{e, e'K^)Y,'^ . These new data sets span a range of momentum transfer 
from 0.7 to 5.4 GeV^, invariant energy W from 1.6 to 2.6 GeV, and the full center-of-mass angular range. 

Our data have been compared to predictions from several available theoretical models that have varying sensitivities 
to the s-channel resonance contributions. The increased statistical precision of these new data will enable improved 
fits either for effective Lagrangian models or for coupled-channels model fits incorporating both photo- and electropro- 
duction data that will be carried out by several groups in the near future 
Center (EBAC) 



62| , including the Excited Baryon Analysis 



63| at Jefferson Laboratory. The analysis of the full set of the world's data in this manner is essential 
to map out the full spectrum of excited states of the nuclcon to better determine the structure of the nuclcon and its 
associated degrees of freedom, both of which are necessary to better understand the strong interaction and QCD. 

The new CLAS A polarization data sets at 4.261 and 5.754 GeV have also been used to extract the longitudinal- 
to-transverse structure function ratio at 6*^™- = 0° in the range from 1.5 to 2.5 GeV^, extending the existing 
CLAS measurements taken at 2.567 GeV near Q^^l.O GeV^. These new data, given the statistical uncertainties, 
could indicate a non-trivial evolution of the structure function ratio in the range from to 2.5 GeV^, that 
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peaks near unity at Q^=1.5 GeV^. These results indicate that extraction of the kaon form factor using the standard 
Chew-Low extrapolation technique can only be carried out in the limited kinematic range where cjl is sizable. 

Finally, the data have been compared to two simple semi-classical partonic models including multi-gluon exchange 
that were designed to account for the strikingly simple phenomenology seen in the kinematic dependence of the 
polarization data. While the two models make very different assumptions regarding the reaction mechanism leading 
to production of the K^K final state and different quantum numbers of the produced ss pair, we have provided 
suggestions for testing them by comparing polarization data for K'^K to K*^ K final states. Disentangling the true 
reaction dynamics in a partonic model is relevant to probe the appropriate quark-pair creation operator that governs 
the transitions to the final state particles and to shed light on the relevance of quark-gluon dynamics in a domain 
thought to be dominated by mcson/baryon degrees of freedom. 

We would like to acknowledge the outstanding efforts of the staff of the Accelerator and the Physics Divisions at 
Jefferson Lab that made this experiment possible. This work was supported in part by the U.S. Department of Energy, 
the National Science Foundation, the Italian Istituto Nazionale di Fisica Nucleare, the French Centre National de 
la Recherche Scientifique, the French Commissariat a I'Energie Atomique, and the Korean Science and Engineering 
Foundation. The Southeastern Universities Research Association (SURA) operated the Thomas Jefferson National 
Accelerator Facility for the United States Department of Energy under contract DE-AC05-84ER40150. 
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